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1 Introduction 

The goal of this paper is to present a topological method, which allows to 
establish, in finite computation, the existence of an infinite number of symmetric 
periodic orbits for dynamical systems with a reversing symmetry. 

The role and impact of time reversing symmetries in dynamical systems has 
been extensively covered in the literature, see |Lal ILall IDll ID2j and references 
given there. In this paper we will restrict ourselves to the discrete time case. 
Before we outline our results we need first to introduce a few definitions. 



^research supported by Polish State Committee for Scientific Research grant 2 P03A 041 

24 



1 



Definition 1 An invertible transformation M : Q — > Q is called a symmetry 
of the discrete dynamical system induced by f : Q — > Q if 

Mo f = f o M. 

It is easy to see that a symmetry of a discrete dynamical system transforms 
trajectories into trajectories. 

Alternatively, a dynamical system may admit transformations S : 51 — > 17 
that map trajectories into other trajectories reversing the time direction. 

Definition 2 JjLalf An invertible transformation S : fl — > is called a re- 
versing symmetry for the discrete dynamical system induced by f : CI — > fl 
if 

Sof = f-^oS, 
or equivalently S^^ o/oS'o/ = Id|o. 

Definition 3 Let (p : T x H. fl (where T = M. or T = Z) be a dynamical 
system. For each point xq G fl, we define its orbit r(a::o) by 

T{xo) — {x \ X — if(t, Xq) for some t G T} 

To describe the symmetry properties of orbits, following jGSSL ILal| we in- 
troduce the notion of isotropy subgroup. 

Definition 4 Consider a dynamical system on a phase space U, with a symmetry 
(reversing symmetry) S. 

If S'(r(a;o)) — r(a;o) for some xq G il, then we say that the orbit r(a;o) is 
S-symmetric. If the symmetric orbit r(a;o) is periodic, then we will say that xq 
is a symmetric periodic point. 

If S is clear from the context, then we will often drop it and speak of a 
symmetric orbit. 

Let P : X X he a map with a reversing symmetry S. Let Fix (S) = 
{x I S{x) — x} be the fixed point set for S. 

A standard method for studying symmetric periodic orbits for systems with 
a reversing symmetry is the Fixed Set Iteration (FSI) method |LQ| ILalj (also 
known as DeVogelaere method |DV| 'l. This method is based on the intersections 
of the iterates of Fix (S) and Fix {PoS) (see [Lall Prop. 1.2.2] for more details). 
Below we give its simplified version. 

Theorem 1 //yi e Fix (5) and P^{yi) 6 Fix (5*) for some k G Z+, then yi is 
a symmetric periodic point for P and its principal period divides 2k. 

Theorem ^ shows that to obtain a symmetric periodic orbit of period k, we 
have to examine the fc/2-iterate of P, which is an apparent obstacle for obtaining 
infinite number of 5-symmetric periodic orbits with unbounded periods in a 
finite computation. 
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One way to overcome this problem is the method proposed by Devaney 
in |D2| . To describe the Devaney 's method, assume that P : R^" R^" 
is a C-'^-map with an reversing symmetry S and Fix (S) is an n-dimensional 
manifold. Now, if we have p G Fix (5), which is a fixed (periodic) point for P 
and the unstable manifold of p, W'^ {p), intersects transversally Fix (5), then by 
symmetry W^{p) = S{W^{p)), hence W^{p) intersects Fix (5) transversally at 
q. Now we apply P'^ to the neighborhood of q and it is easy to see, that for 
k large enough there exists a disk Dk in Fix (5*), such that P^{Dk) intersects 
Fix {S) transversally in the vicinity of p. Now by Theorcmnoi^^ obtains periodic 
points of arbitrary high periods. 

The method proposed in this paper bears some resemblance with the De- 
vaney's method, as it also is based on some kind of transversality, which is 
preserved under the iteration of the map. Our method is purely topological and 
is based on the notion of covering relation (see |GiZ| 'l as a tool for the prop- 
agation of the topological transversality. Here we will informally outline our 

p 

approach. We say that a cube N P-covers a cube M, denoted by N =^ AI, if 
the image of the cube N under the map P is stretched across the cube M in a 
topologically nontrivial manner (see Definition ^ . These cubes together with 
the corresponding choices of coordinate systems will be referred to as h-sets (the 
letter h suggesting the hyperbolic- like directions) . Now in each h-set N we can 
define horizontal and vertical disks (see Definitions IIUI and II 1|) . In Section |3| we 
prove (see Theorem that if we have the chain of covering relations 

N = No ^ Ni ^ . . . ^ Nk ^ M, (1) 

then for any horizontal disk H in Nq and any vertical disk V in Nk , there exists 
X e H, such that P*(a;) e TV; for i = 1, . . . , fc and P''{x) e V. Now if Fix (S) 
forms a horizontal disk in N and Fix (S) forms a vertical disk in M, then from 
Theorem^ we obtain symmetric periodic point. Observe that if we can build 
a chain of covering relations ^ linking N and M of arbitrary length, then we 
will have symmetric periodic points of arbitrary high period. For example it 
was shown that this happens for suitable 2-dimensional Poincare maps for the 
Michelson system arising from the Kuramoto-Sivashinsky PDE |W|, the planar 
restricted three body problem modelling the motion of the Oterma comet in the 
Sun- Jupiter system ,WZ2) or the Henon-Heiles Hamiltonian |AZ) . In the above 
mentioned applications we had one unstable and one stable directions. In this 
paper we apply our method to a four-dimensional reversible map to prove the 
existence of symbolic dynamics and an infinite number of symmetric periodic 
orbits. The main feature, which makes this example interesting, is the fact that 
both stable and unstable directions are two-dimensional. The proof is computer 
assisted, i.e., rigorous numerics is was used to verify assumptions of abstract 
theorems. 

The proposed method was introduced in W in the planar case and direct 
coverings. In this case the proof of the transversality theorem (corresponding 
to Theorem 13 is very simple and is based on the connectivity argument. Un- 
fortunately this proof cannot be generalized to higher dimension or to include 
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coverings induced by inverse mappings. 

The content of this paper can be described as follows. In Section|21we define 
topological notions: h-sets, covering relations and backcoverings relations. In 
Section 13 we prove the main transversality theorem (Theorem IS} . In Section 0] 
describe how it can be applied to reversible dynamical systems in general. In 
Section |31 we consider an application of our method to a four-dimensional map 
with an reversing symmetry. In Section|Blwe describe how to verify the existence 
of covering relations by computer. 

2 Topological tools: h-sets and covering rela- 
tions 

In this section we present main topological tools used in this paper. The crucial 
notion is that of covering relation [GiZj. 

2.1 h-sets 

Notation: For a given norm in M" by Bn{c, r) we will denote an open ball of 
radius r centered at c S M". When the dimension n is obvious from the context 
we will drop the subscript n. Let S^^{c,r) = 9i3„+i(c, r), by the symbol 5*" we 
will denote ^"(0, 1). We set MP_^ {0}, Bo(0,r) = {0}, dBa{0,r) = 0. 

For a given set Z , by int Z, Z, dZ we denote the interior, the closure and the 
boundary of Z, respectively. For the map h : [0, 1] x Z ^ M" we set ht — h{t, ■). 
By Id we denote the identity map. For a map /, by dom(/) we will denote the 
domain of /. Let f : fl C M" — + M" be a continuous map, then we will say that 
X C dom(/^^) if the map /"^ : X M" is well defined and continuous. For 
N C fl, iV-open and c G ]R" by deg(/, N, c) we denote the local Brouwer degree. 
For the properties of this notion we refer the reader to ^ (see also Appendix 
in [niZl)- 

Definition 5 JjGiA Definition 1] A h-set, N, is the object consisting of the 
following data 

• \N\ - a compact subset o/R" 

• u{N), s{N) e {0, 1,2,...}, such that u{N) + s{N) = n 

• a homeomorphism cn : K" ^ K" = K"^^^ x such that 

CNm) = B^){^. 1) X S:^(0, 1). 

We set 

7Vc = b;:^(o,i) X 51^(0,1), 
N- = 5^:^(0, 1) X s:^(o, 1) 

7V+ = b:^)[o, 1) X as:^(o, 1) 

N-^c-^\N~), N+^c-^\N+) 
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Hence a /i-set, N, is a product of two closed balls in some coordinate system. 
The numbers, u{N) and s{N), stand for the dimensions of nominally unstable 
and stable directions, respectively. The subscript c refers to the new coordinates 
given by homeomorphism cjv- Observe that if u{N) = 0, then N~ = and if 
s{N) = 0, then TV"*" = 0. In the sequel to make notation less cumbersome we 
will drop the bars in the symbol |A^| and we will use TV to denote both the h-sets 
and its support. 

Definition 6 IGiZi Definition 3] Let N he a h-set. We define a h-set N'^ as 
follows 

• l^^l = 1^1 

• u{N'^) = s{N), s{N^) = u{N) 

• We define a homeomorphism c^t : M" ^ R" = R"(^^) x W<-'^'^\ by 

cnt{x) = j(cjv(a;)), 
where j : R''^^) x W'-^'^ R"'^) x M"(^) is given by j{p, q) = {q,p). 



Observe that N'^'^ = N and N'^' — N^. This operation is useful in the 
context of inverse maps. 



2.2 Covering relations 

Definition 7 lGi2X Definition 6] Assume that N, M are h-sets, such that u{N) — 
u{M) — u and s{N) = s{M) = s. Let f : N M" be a continuous map. Let 
fc — Cm ° / ° c]^^ : Nc R" x . Let w be a nonzero integer. We say that 

(N f -covers M with degree w) iff the following conditions are satisfied 

1. there exists a continuous homotopy h : [0, 1] x Nc — > R" x R^, such that the 

following conditions hold true 

ho = fc, (2) 

h{[0,l],N-)nMc = 0, (3) 

h{[0,l],Nc)nM+ = 0. (4) 

2. There exists a map A : M" M", such that 

hi{p,q) = (A(p),0), /or 5^(0,1) anrfge 57(0,1), (5) 
A(aS„(0,l)) C R"\ 57(0,1). (6) 

Moreover, we require that 

deg(A,S7(0,l),0) = «;, 
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Note that in the case u ^ 0, if N M, then /(|iV|) C int |A/| and w = l. 

Intuitively, N M if / stretches N in the 'nominaUy unstable' direction, 
so that its projection onto 'unstable' direction in M covers in topologically 
nontrivial manner projection of M . In the 'nominally stable' direction N is 
contracted by /. As a result N is mapped across M in the unstable direction, 
without touching M+. It is also very helpful to note that the degree w in the 
covering relation depends only on ^|aB„(o,i)- 

Definition 8 \Gi2\ Definition 7] Assume N, M are h-sets, such that u{N) = 
u{M) = u and s{N) = s{M) = s. Let g : Q C M" K". Assume that 
: \M\ M" is well defined and continuous. We say that N ^= M (N 

g-backcovers M with degree w) iff M'^ 7V^. 
The following theorem was proved in |GiZj . 

Theorem 2 \GiZ[ Theorem 9] Assume Ni, i = 0, . . . ,k, N/^ = Nq are h-sets 
and for each i = 1, . . . , k we have either 

7V._i (7) 

or \Ni\ C dom(/-^"^) and 

7V,_i N,. (8) 
Then there exists a point x G int \Nq\, such that 

/,o/,_io-- -0/1(2;) e mt\N^\, i^l,...,k (9) 
/fc o /fc-i o • • • o /i(a;) = X (10) 

Obviously we cannot make any claim about the uniqueness of x in Theo- 
rem El 

Theorem 121 shows that both the direct and the inverse covering relation can 
be treated on the same footing. This justifies the following definition. 

Definition 9 Assume iV, M are h-sets and P is a continuous map. We say 
that 

If one of the two following conditions is satisfied 

N C dom (P) and N =^ M 
\M\ C dom [P-^) and N ^ M. 

P w 

We would like to stress, that the relation N <;==4> M is not symmetric. 
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3 The topological transversality theorem 



The goal of this section is to state and prove the main topological transversality 
theorem for chain of covering relations. For this end we need first to define the 
notions of vertical and horizontal disks in an h-set. 

Definition 10 Let N be an h-set. Let b : i?„(7v)(0, 1) — > |iV| be continuous 
and let be — o b. We say that b is a horizontal disk in N if there exists a 
continuous homotopy h : [0, 1] x _B„(7v)(0, 1) Nc, such that 

ho = foe (11) 

hi{x) = ix,0), for allxeB;^j{0,l) (12) 

h{t,x) e N~, for all t e [0,1] and X e dR^{0,l) (13) 

Definition 11 Let N be an h-set. Let b : Bs(n){0, 1) \N\ be continuous and 
let be = cn o b. We say that b is a vertical disk in N if there exists a continuous 
homotopy h : [0, 1] x i?s(Ar)(0, 1) Nc, such that 

ho = be 

hi{x) — (0,2;), for all X E Bg^ff^ (0,1) 
h{t, x) € N+, for all t G [0, 1] and x G dB^){0, 1). 

It is easy to see that b is the horizontal disk in iff is the vertical disk in iV^. 

We would like to remark here that the horizontal disk in TV can be at the 
same time also vertical in N. An example of such disk is shown on Fig. ^ 
In case homotopies used in the definitions of horizontal and vertical disks are 
different. The existence of such disks, which are both vertical and horizonal 
will play very important role in our method for detection of an infinite number 
symmetric periodic orbits for maps with reversal symmetry. 




Figure 1: The curve b is both horizontal and vertical disk in N. In this example 
u{N) = s{N) = 1. 
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or 



Now we are ready to state and prove the main topological transversality 
theorem. A simplified version of this theorem was given in W, for the case of 
one unstable direction and covering relations chain without backcoverings. The 
argument in |W], which was quite simple and was based on the connectivity 
only, cannot be carried over to larger number of unstable directions or to the 
situation when both covering and backcovering relations are present. 

Theorem 3 Let k > I. Assume Ni, i = 0,...,k, are h-sets and for each 
i = 1, . . . ,k we have either 

iV,_i ^m>' (14) 

iV,_i N,. (15) 

Assume that bo is a horizontal disk in Nq and be is a vertical disk in Nk- 
Then there exists a point x G int |A'^o|; such that 

X = bo(t), for some t € Bu(No)iO,^) (16) 
/.o/,_io...o/i(x) G intlA^.I, i = l,...,k (17) 
fko fk-i° ■ ■ ■ ° h{x) = be{z), for some z e Bs^Nt){OA) (18) 
Proof: Without any loss generality we can assume that 
CNi = Id, for i = 0, . . . , k. 

Then 

fi^fi.c, for « = 1, . . . ,fc, 
N, = iV,,„ 7V± = iV± for i = 0, . . . , A: 

We define gi — fj^^, for those i for which we have the back-covering relation 

Notice that from the definition of covering relation, it follows immediately 
that there are m > 0, s > 0, such that u{Ni) — u and s{Ni) — s, for all 
i = 0, . . . ,k. 

The idea of the proof is to rewrite our problem as a zero finding problem for 
a suitable map, then to compute its local Brouwer degree to infer the existence 
of a solution. 

As a tool for keeping track of the occurrences of coverings and backcoverings, 

we define the map dir : {1, . . . , fc} {0, 1} by dir («) = 1 if Ni^i ^^=^ Ni and 

dir {i) = if Ni-i ^4=^ Ni. For « = 1, . . . , /c let /li be a homotopy map from the 

definition of covering relation for A^^-i Ni or A^^-i <J=' Ni. In the case of 
a direct covering (i.e. dir (i) = 1), the homotopy hi satisfies 

h,{0,x) = fi{x), wherea;GM"+^ (19) 

h,{l,{p,q)) = (A,(p),0), where p G M" and g G R^ (20) 

/i,([0,l],7Vr JniV, = 0, (21) 

h,{[0,l],N,^i)nN+ - 0. (22) 
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In the case of a backcovering (i.e. dir (i) = 0), the homotopy hi satisfies 

h^{0,x) = g,{x), where a; e R"+^ (23) 

h,{l,{p,q)) = {0,A,{q)), where p e M" and g G R", (24) 

/ii([0,l],iV+)nA^,_i = 0, (25) 

h,{[0,llN,)nNr_^ = 0. (26) 

Let ht and hz be the homotopies appearing in the definition of the horizontal 
and vertical disk for fep b^, respectively. 

It is enough to prove that there exists t G 3^(0,1), z e -6^(0, 1) and Xi G 
int Ni for i = 1, . . . , /c such that 

6o(i) = a;o, 
fi{xi-i) = Si, ifdir(i) = l, 

gi{xi) = if dir (i) = 0, (27) 

be{z) = 

We will treat (|27|l as a multidimensional system of equations to be solved. 
To this end, let us define 

n = B;:(0, 1) X A^i X • • • X Nk-i X 57(0, 1) 

A point a; G n will be represented by a; = {t,xi,. . . ,Xk~i,z). 

We define a map F = (Fi,...,Ffe) : H ^ Ri-^+^)'' as follows: for i = 
2, . . . , fc — 1 we set 



Fi{t,xi,...,Xk-i,z) 
For z = 1 we set 

Fi(t,xi,. . . , Xk-i, z) = 
For i = k we define 

Fk{t,xi, . . . ,Xk-i,z) = < 



x^ - fi{xi-i) ifdir(i) = l, 
a^j-i - Oiixi) if dir {i) = 0. 



Xi~fi{bo{t)) ifdir(l) = l, 
bo{t)-gi{xi) ifdir(l) = 0. 



beiz) - fk{xk-i) ifdir(fc) = l, 
Xk-i - gk{be{z)) ifdir(fc)=0. 



With this notation, solving the system (|27|l is equivalent to solving the equa- 
tion F{x) = in intn. 

We define a homotopy = (ffi, . . . , iJfe) : [0, 1] x H ^ as follows. 

For i = 2, . . . , fc — 1 we set 



Hi{X,t,xi, . . .,Xk-~i,z) 



Xi - hi{X,Xi^i) ifdir(i) = l, 
Xi-i — hi{X, Xi) if dir {i) = 0. 
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For z = 1 we set 

Hi{X,t,xi,...,Xk-i,z) 
For i = k we define 

HkiX,t,xi, . . .,Xk-i,z) = 



xi - hi{\ht(\,t)) ifdir(l) = l, 
ht{\t)-hi{\,xi) ifdir(l) = 0. 



/i2(A, z) - /ife(A, a;fe_i) ifdir(fc) = l, 
Xk-i - hk{X,h^_{X, z)) ifdir(fc)=0. 



Notice that H{0, x) — F{x). The assertion of the theorem is a consequence 
of the following two lemmas, which will be proved after we complete the current 
proof. 

Lemma 4 For all A e [0,1] the local Brouwer degree deg(_ff(A, •), int H, 0) is 
well defined and does not depend on A. Namely, for all A G [0, 1] we have 

deg{H{X, •), int H, 0) = deg(i/ (1, •), int H, 0). 

Lemma 5 

|deg(H(l,-),intn,0)| ^\wi-W2 Wk\ 

We continue the proof of Theorem O Since F = H{Q,-), from the above 
lemmas it follows immediately that 

deg(F, int n, 0) = deg(iJ(0, •), int H, 0) = deg(iJ(l, •), int H, 0) 7^ 0. 

Hence there exists a; € 11 such that F{x) =0. 1 

Proof of Lemma |4l From the homotopy property of the local Brouwer degree 
(see Appendix in jGiZp it is enough to prove that 

H{X, x) ^ 0, for aU x e d{U) and A £ [0, 1]. (28) 

In order to prove (|28|l . let us fLx x — {t, xi, . . . , Xk-i, k) G 911. It is easy to 
see that one of the following conditions must be satisfied 

t e 9B„(0,1), (29) 

z e dBs{0,l), (30) 

Xi e N^, for some i = 1, . . . , A; — 1, (31) 

Xi e N~ for some i = I, . . . ,k — I. (32) 

We will deal with all above cases separately. 

Consider first (EH- Let us fix A e [0,1]. Let xq = ht{X,t). From Def. [TUl 
it follows that xq G . There are now two possibilities: either dir (1) ~ 1 
(direct covering) or dir(l) = (backcovering) . Assume that dir(l) = 1. From 
condition (|^ it follows that /ii(A,xo) ^ iVi, hence Hi{X,x) ^ 0. Assume now 
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that dir (1) = 0. We have ht{X, t) 6 and since by ^ hi{X, Ni) n = 0, 
hence i?i(A, a:) ^ 0. 

Consider now H3U|) . Let us fix A G [0, 1] and let Xk — hz{X, z). From Def.llllit 
follows that Xk & . Now if dir (k) — 1, then from condition (|22|l /ife(A, A^fc_i)n 
Ar+ = 0, hence HkiX,x) ^ 0. If dir(fc) = 0, then from ^ if follows, that 
hk{X,Xk) i- Nk-i, hence Hk{X,x) ^ 0. 

For each of cases (|31|l and (|32() we have to consider the following four pos- 
sibilities 







(33) 






(34) 






(35) 




4^ n/4^n^+i. 


(36) 



Assume first that Xi £ Nl^^ . If H33|l or (|34() holds true, then from H22() we obtain 

hi{X, Xi-i) 7^ Xi, 

for every A £ [0, 1] and every Xi^i € A^i-i- If or is satisfied, then from 
(123 it results that 

hi{X, Xi) ^ Xi-i, 

for every A G [0,1] and every Xi^i e jA^i-ij. This proves if Xi £ A^^^, then 
H{X,x) ^ for any A G [0,1]. 

Assume now that Xi £ 'N^ . If or holds true, then from ||2Ill it 
follows that for every A G [0, 1] and every Xi^\ G NiJ^\ we have 

K^\{X,Xi) ^ Zi+i. 

If (|34|l or H36|) is satisfied, then from H26() we obtain 

/ii+i(A, 7^ Xi, 

for every A G [0, 1] and every x^+i G A^i+i- This proves that if Xi G N~ , then 
F(A,x) 7^ for any A G [0,1]. I 

Proof of Lemma [Sj Let us represent Xi for i — 1,...,A; — 1 as a pair 
— iPiiQi)i where pi G and G M''. In this representation the map 
iJ(l,t,Pi,gi, . . . ,pk-i,qk^i,z) . . . ,Pk,qk) has the following form (for 

a = 0) 

• if i = 2, . . . , A: — 1 then 

if dir (i) = 1, then = (1 - - Ai(pi_i), Qi = qt, (37) 

if dir (i) = 0, then pj=pj_i, = (1 - Q;)gj_i - yli(5i) (38) 
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• if i = 1, then 

if dir (1) ^ 1, then pi (1 - a)-pi - Ax{t), qi = qi (39) 
if dir (1) = 0, then pi = t, qi = -Ai{qi) (40) 

• li i = k, then 

if dir (fc) = 1, then pk = -Ak(pk-i), qk ^ z (41) 
if dir {k) = 0, then pk = Pk^i, qk = (1 - a.)qk-i - Ak{z) (42) 

The above equations define a homotopy C : [0, 1] x 11 ^ We wiU show 

that deg(C(Q;, •),intn, 0) is independent of a and then we compute the degree 
ofC(l,-). 

Lemma 6 For any a G [0, 1] 

deg(C(a,-),intn,0) = deg(C(l, •), int H, 0). 

Proof: From the homotopy property of the local degree (see Appendix in |GiZ| l. 
it follows that it is enough to prove that 

C(a, x) ^ 0, for aU x £ an and a G [0, 1]. (43) 

Let us take x = (i,pi, gi, . . . ,pk-i, qu-i, z) G 911. One of the following condi- 
tions holds true 



t 


G 


5", 








z 


G 


^^ 








Pi 


G 




for some i — 


2,. 


..,k-l 


qi 


G 


^^ 


for some i = 


2,.. 


.,k-l. 



Assume that t G S"". If dir (1) = 1, then \\Ai{t)\\ > 1, hence ||pi|| > \\Ai{t)\\ ~ 
IIpiII > 0. If dir (1) = 0, then pi = t ^ 0. 

Assume that z G S''. If dir (k) — 1, then = z 7^ 0. If dir (fc) = 0, then 
\\Ak{z)\\ > 1 and we obtain ||gfe|| > \\Ak{z)\\ - \\qk-i\\ > 0. 

Assume that pi G S"". If dir (i + 1) = 1, then pi+i ^ 0, because from 
condition ((HJ it follows that 

||Am(ft)ll>l>ll(l-«)mi||, (44) 

for any pi+i G B„(0, 1). 

If dir (i + 1) = 0, then obviously pi+i = pi ^ 0. 

The argument for the case qi G S"" is similar. 1 

Now we turn to the computation of the degree of C(l,-). Observe that 
C(l, •) has the following form: for i = 2, . . . ,k — 1 

Pi = ~Ai{pi_i), qi = qi if dir (i) = 1 (45) 
Pi = Pi-i, qi = -A.,{q,) if dir (i) = 0, (46) 
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for i — I 



and for i = k 



Pi^-Ai{t), qi^qi ifdir(l) = l (47) 
Pi -= t, qi = if dir (1) = 0, (48) 



Pk = -Ak{pk-i), qk^z if dir (1) = 1 (49) 
Pk^t, qk = -Ak{z) ifdir(l) = 0. (50) 

From the product property of the degree (see Appendix in |GiZ| ') it foUows 
that 

_ |deg(C(V),n,0)| = 

|n,;6dir-i(i) deg(-A,, b;:(0, 1), 0) • n,gdir-i(o) deg(-A„ Sl(0, 1), 0)| . 

In the formula above if dir (j) = (for j = 0, 1), then the corresponding prod- 
uct is set to be equal to 1. Similarly if u = or s = 0, then the corresponding 
product is also set equal to 1. 

From CoUorary 18 in |GiZ| it follows that 



deg(-A, U, 0) = (-1)" deg(A, [/, 0). (51) 
This finishes the proof. I 

4 Reversing symmetry and covering relations. 

In this section we apply the tools developed in previous sections to the study of 
symmetric periodic orbits for reversible maps. 

Theorem 7 Let S be a reversing symmetry for the local dynamical system in- 
duced by the map P. Assume that 

Mo ^ A/i ^ M2 • • • ^ M„_i 4U M„, 

and 

1. there exists a horizontal disk in Mq contained in Fix (S*) 

2. there exists a vertical disk in M„ contained in Fix (S) 

Then there exists x G Mq, such that 

S{x)^x, P2"(x) 
the orbit of x is S -symmetric 
P\x) ^ Mi for i = 1, ... , n, 
P"-+\x) G S'(M„-,;) /or i = 1, . . . , n. 
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Proof: From Theorem O it follows that there exists x £ Fix (S*), such that 

P\x) e Mi for j = 1, . . . , n, 
P"{x) e Fix(S'). 

The assertion now follows from the reversing symmetry of P. i 

Now we turn our attention to the action of symmetry on h-sets and covering 
relations. 

Definition 12 Let N be an h-set mR". Let L : R" be a homeomorphism. 

We define an h-set L * N as follows 

• |L*7V| = L{\N\), 

• u{L ^N)^ u{N) and s{L x= A^) = s{N), 

• Cl*N = Cat o L^i. 

We define an h-set L'^ * N by 

L^ *N = {L^Nf 

Informally speaking, L* N is just a natural symmetric image of N and L'^ * N 
is the symmetric image of N, but we additionally switch the 'expanding' and 
'contracting' directions. 
We have the following 

Lemma 8 Let S : R^" — > R^" be a reversing symmetry for a map P and 
N <^ M, then 5^ * M 4^ S'^ *N. 

Proof: From the definition of covering relations and reversing symmetry it 
follows immediately, that 

if N M, then * M <^ * N 
if N <^ M, then 5*^ * M S'^ * N 

I 

The reversing symmetry maps, in a natural way, horizontal disks to vertical 
disks and vice versa. Namely, we have the following obvious lemma. 

Lemma 9 Let S : R" — > R" a homeomorphism, N be an h-set and j be a 
horizontal (vertical) disk in N. 

Then S{j) is a vertical (horizontal) disk in S'^ * N . 

In the context of proving the existence of an infinite number of symmetric 
periodic orbits symmetric h-sets are of special importance. 

Definition 13 Let S be a reversing symmetry. We say that an h-set N is 
^^-symmetric if S"^ * N = N . 
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It is easy to see, that if N is 5"^-syinnietric h-set, then u{N) — s(N) and the 
dimension of the phase space n = u{N) + s(N) must be even. 

The following theorem, which is an easy consequence of Theorem [7\ illus- 
trates our method of proving of the existence of an infinite number of symmetric 
periodic orbits. 

Theorem 10 Let S be a reversing symmetry for the local dynamical system 
induced by the map P. Let Mi for z = 0, 1, 2, 3 be h-sets and Mi n Mj = for 
i ^ j ■ Assume that 

Mo 4^ Ml 4^ M2, 

Mo ^ Mo 
Mq and M2 are S"^ -symmetric 
M3 = S*^ * Ml 
there exists a horizontal disk in Mq contained in Fix (S) 
there exists a horizontal disk in AI2 contained in Fix [S) 

Then for any sequence (ap, ai, . . . , a„) £ {0, 1, 2, 3}"+^, satisfying the fol- 
lowing conditions 

aoG{0,2}, a„e{0,2} 
if Ui = 0, then a^+i G {0,1} 
if OLi — 1, then a^+i = 2 
if Ui — 2, then a^+i — 3 
if ai — 3, then a^+i = 1, 

there exists x G lAfctgl, such that 

S{x)^x, 
the orbit of x is S -symmetric 
P\x) e \M,\ fori^l,...,n, 
P"+'(.t) e Si\Mn-^\) fori^l,...,n. 

Proof: From Lemma |51 it follows that in Mq and M2 there exist vertical disks 
contained in Fix(S'). We can now apply Theorem [71 to the chain of covering 
relations 

Af„„ ^ 4^ ... 4^ A/„„ 



5 One four-dimensional reversible example. 

In this section we present the application of the method introduced throughout 
the paper to a four-dimensional reversible map. As a consequence we obtain 
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the existence of chaotic dynamics and the existence of an infinite number of 
symmetric periodic orbits for a certain iteration of such map. The proof is com- 
puter assisted, i.e., rigorous numerics is used to verify assumptions of abstract 
theorems. The main feature, which makes this example interesting is the fact 
that both stable and unstable directions are two-dimensional and the map itself 
is not close to a product of two two-dimensional maps, with unstable dimension 
each. 



5.1 An example of four-dimensional reversible map. 

Let / : R" R", n > be a continuous map, F : R^" — > R^" be a map defined 

by 



F 


X 








y 




X 



and let S{x, y) — (y, x). It is straightforward to show that S o F o S o F — Id. 
Therefore F is a reversible homeomorphism of R^" . In suitable coordinates we 
may rewrite F as follows 



F 



X 




y. 





i2- 



y) 



(52) 



and the reversing symmetry S will be given by S{x,y) = 
Let us fix n = 2 and let / : R^ ^ R^ be defined by 



/ 



,y)- 



Xi 






- Xx) - 


1-4 


- X2 






.0:2(1 


— X2) - 


1-4 


+ Xl 



(53) 



In the remainder of this section we will investigate the map F given by l|52|l . 
where / is as above. The map S{xi,X2, 2/1,2/2) = (— a;i, —X2, yi, 2/2) is a reversing 
symmetry of F. 

It is easy to verify that each solution of 



F{xi,X2,yi,y2) = (2^1, 2:2, 2/1, 2/2) 



(54) 



satisfies xi — X2 — 0, hence all fixed points of F are symmetric. Solving of 
Eq. (|54|l leads to the following system of equations 

[2/? + (2/2 + 1)2 = 9 
[(2/1 + 1)2-2/2 = 1 

describing an intersection of a circle with a hyperbola in four points as was 
shown in Fig. 12 Fixed points Pi,P2 are hyperbolic with real eigenvalues, the 
fixed point P3 is hyperbolic with four complex eigenvalues, the point P4 possesses 
two complex eigenvalues on the unit circle and two real eigenvalues, i.e., it is of 
elliptic-hyperbolic type. 
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Figure 2: The location of the fixed points of F. 

The fixed points may be exactly computed (for example using Mathematica). 
However, for our further consideration it is sufficient to use two approximate 
fixed points, which we will still denote by Pi, given by 

Pi = (0,0,-2.9288690017630725,-1.649404627725545), (55) 
Pa = (0,0,2.199939462565084,-3.0396731015162355). (56) 

We will show that in the vicinity of Pi and P2 the map F'^ has symbolic dynamics 
on two symbols and there exist an infinite number of symmetric periodic points 
of an arbitrary large period. 

Before we proceed with the statement of main results for F we need to 
discuss how we represent h-sets. 

5.2 Representation of h-sets in M". 

To define a h-set N we need to specify a homeomorphism cjv of M" and two 
numbers u{N) and s{N) (see Def. IS)). Since we will use the computer in order 
to verify covering relations, the homeomorphism must be representable by the 
machine. The simplest case is to take an affine map. We will use the maximum 
norm on M", i.e. ||a;|| = max^ \xi\ and we treat vectors as columns. 

Let x,ui,..., Uk, si, . . . , Sn-k G R", < A: < n be such that the vectors 
ui, . . . , Ufc, si, . . . , Sn-k are linearly independent. We define a matrix M £ R"^" 

by 

M = [Ui, . . . Si, . . . , Sn-k]- 

We define a h-set 

iV = l){x,Ui, . . . ,Uk,Si, . . . ,Sn-k) 
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as follows 

u{N) = k, 

s{N) = n-k, 

\N\ = M(S;(0,l)) + x = M(S;:^(0,l) x;B:^(0,l)) + a:, 

cn{v) = M-'^{v-x), forweM" 

Hence, the h-set N defined above is a parallelepiped centered at x. 

In the sequel we will work in K^*^ with s{N) — u{N) = k. In this case we 
will use also the notation N = [){x,M), where M is a linear isomorphism of 
R^*^ . The first k columns of M correspond to unstable directions and the last k 
columns of M correspond to stable directions. 



5.3 Important h-sets and covering relations between them 



As was mentioned before Pi and P2 are good numerical approximations to two 
hyperbolic fixed points with two-dimensional stable and unstable manifolds. We 
choose Uj,Sj S K^, J, j = 1, 2 to be a good numerical approximations of unstable 
and stable eigenvectors of DF{Pi). Put 

^0.527847408170044 
0.254065286036574 
0.730261232439584 
0.351491787265563 

"0.233876807615845 
0.485903716548415 
0.365235930520818 
^0.758816138574061 

J — 



-0.05726452423754 ' 
0.594572575636284 
-0.0768865282444865 
0.7983061369797889 

0.8918103319483236 
-0.0858921121857865 

0.4421352370808943 
-0.0425829663821858) 



(57) 



and put si — S[ul), i,j — 1, 2. We define matrices Mi for i = 1, 2 by 

Ml = 0.012[ul,ul,sl,sl], M2 = 0.31[ul,ul,sl,sl]. 
We define two S'^ symmetric h-sets centered at Pi and P2 by 

iVi =f)(Pi,Afi), N2 = HP2,M2). 



(58) 



(59) 



The important remark is that the sets \Ni \ and |A^2| are disjoint. The projection 
of |A^i| and \N2\ onto (2/1,1/2) coordinates is presented in Fig.O The following 
lemma was proved with a computer assistance. 

Lemma 11 The following covering relations hold: 



Ni ^ Ni, 



N2 ^ N2. 



The details of the proof will be presented in Sectional 
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-3-2-10 1 2 



Figure 3: The projection of the sets \Ni\ and \N2\ onto (1/1,2/2) coordinates. 

Remark 12 From above lemma and Theorem 2 it follows immediately, that 
there exists fixed points in Ni and N2 ■ However, we cannot claim that they are 
unique in Ni nor S-symmetric. 

Now wc will construct a dynamical link between sets A^i and A^2- More 
precisely, we construct a chain of covering relations connecting Ni and N2. For 
this purpose we look for a point Qi in the neighborhood of Pi on the unstable 
manifold of Pi and such that F'^(Qi) (here k = 10) is close to P2. Then we 
define additional points Qi by taking some forward iterates of Qi. 

To be specific we set 

Qi = 0.0330092432m} - 0.048949m^ + 0.0004931s}, 

Q2 = F\Qi), 
Qs = F{Q2). 

We have 

F-i(Oi) e |iVi|, 

||J^-i(Qi) - Pill < 0.006, 
||i=^'°(Qi)-P2|| < 0.001. 

The points Qi, Q2 and Qa will be used as the centers of new h-sets. We define 

Hi = l)(gi,.001u},.00175u4-005s}, 0.0054), 
H2 = l)(g2,. 28^2,. 28^2, .2s2,.38s2), 
= i}{Q:i,.lhul,.lbul,.l2slA2sl). 



Lemma 13 The following covering relations hold 

F 1 F*— 1 F— 1 F— 1 

Ni ^ Hi ^ H2 ^ H3 ^ N2. 
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The details of the proof will be presented in Section [HI 

Let us comment briefly about the spatial relations between Ni and Hj . The 
set Hi is close to Ni, this is the reason for using the same stable and unstable 
directions in the definitions of these sets. Sets H2 and H3 are close to A^2 and as 

in the previous case we used the same stable and unstable directions for them. 

pi -1 

The numerical evidence of the existence of covering relation Hi H2 is 

presented in Fig. 0] 



4 
2 

-2 
-4 



-7.5 -5 -2.5 2.5 5 7.5 




Figure 4: The numerical evidence of the existence of covering relation Hi 
H2, (left) the set F^{{Hi)'^) projected onto unstable directions of (i?2)c is 
outside the miit ball, (right) the set F^{{Hi)c) projected onto stable directions 
of {H2)c is inside the unit ball. 



5.4 Chaotic dynamics of F. 

Theorem 14 The discrete dynamical system induced by the map IS semicon- 
jugated with the full shift on two symbols, i.e. for an arbitrary {ij)ji=.z 6 {1, 2}^ 
there exists a point xq G l-^iol such that 

f7^(xo) e |7V,j, jez. (60) 

Moreover, if the sequence {ij)j^z is periodic, then the point may be chosen as a 
periodic point for F"^ with the same principal period. 

Proof: From Lemma [TSl we obtain that 

, ^ F. 1 , , . — 1 , , F. — 1 , , F. — 

Ni^ Hi =^ H2 ^ H3 =^ A^2- (61) 

Since the h-sets A'^i and N2 are symmetric (by their definition), the reversing 
symmetry property of F implies 

7V2 = 5^ ★ N2 ^ * H-s ^ * H2 
5*^ ★ H2 V= S'^ i.Hi -U^ S'^ irNi^ Ni. 
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From Lemma fm we get 

7 times 



N,^N,^---^m, (63) 

7 times 



N2-^ N2 =^ ■■■ =^ N2. (64) 

Let (jj)jGZ G {1, 2}^ be a periodic sequence od symbols, i.e. ij+fc = ij for j G Z 
and a certain /c > 0. Let 

^N,,---^N,,^ TV,, (65) 

be a periodic sequence of covering relations, where by Ni. we mean 

a corresponding sequence H61I64|I . Now, Theorem |2 implies that there exists a 
fc-periodic point xq G \Nig \ for F"^ , such that assertion (|60() is satisfied. 

Let {ij)jei G {1:2}^ be a nonperiodic sequence of symbols. We define the 
periodic sequences 

(«j )jez = {•••, {i-fc, • ■ ■ , «o, ■ • ■ , */c}, {i-fe, • ■ • , »o, • ■ ■ , U-}, • ■ ■} e {1, 2}^ 

where — ioi k > 0. Now, for fc > we can find (2fc + 1)— periodic point Xk 
for i^^ such that F'^^ (xk) € l-^i^l, fc > 0, j G Z. Since Nig is a compact set, 
we can find a condensation point xq G |-/Vi„| of {xk}k>o- Obviously, xq satisfies 
assertion (|Sn)l. I 



5.5 Symmetric periodic points for F. 

Theorem 15 There exists an infinite number of symmetric periodic points for 
F with an arbitrary large principal periods. 

The proof of Theorem 1151 is a direct consequence of the following lemma. 
Lemma 16 Let 

be a sequence of covering relations, where 

Vo,Vk G {Ni,N2}, 
V, G {iVi, Hi,H2, H3}, /or J = 1, . . . , fc - 1, 
gi ^ F or gi ^ F^ , for i = 0, . . . , k - I. 

Then there exists a symmetric periodic point xq of F , such that 

xo G |K)|nFix(5), 
ifji o ■■ - ogi o go){xo) G i = 0, . . . , fc - 2, 

(g/c-l o • ■ • 051 o gQ){xo) G \Vk\ nFix(S') 
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Proof: Recall, that the sets Ni , N2 are symmetric and defined by vectors ufjsf, 
i, j = 1, 2 - see(|S7|l and For i = 1, 2 we define a map bi : 'B^{0, 1) — > \Ni\ 

by 

bi{p,q) = AU[p,q,p,q]'^ + Pi, i = l,2. 

We will show, that hi is a horizontal disk in Ni. 
First observe that 

bi,c = CN, o b^{p, q) = {p, q,p, q). (66) 

We define the homotopy hi by 

ht{t,p,q) = (1 - t)p, (1 - t)q). 

It is easy to see, that conditions from DefinitionllOlare satisfied. Namely, 

we have 

hz{0,p,q) = {p,q,p,q) ^ {b^)c{p,q), for (p, gr) £ ^2(0, 1), 
h,{l,p,q) = (p,g,0,0),_for (p,g)^B^(0,l), _ 

e (iVO,~ -95^(0,1) X 5^(0,1), for (p,<7) £9:8^(0, [0,1]. 

This proves that bi is a horizontal disk in Ni. 

Let us remind the reader that S{Pi) = Pi and S'(up = s^. Hence we obtain 

S{br{p, q)) = S{h{pu\ + qu^ + ps\ + qs\) + P,) = 
= h {pS{u\) + qS{u\) + pS{s\) + qS{s^)) + S{P,) - b,{p, q). 

where fci = 0.012, k2 = 0.31 are the coefficients used in (|58f) to define N^. This 
proves that bi are contained in Fix (5*) . 

We will show that bi is also vertical disk in Ni. Namely, if follows from 
Lemma El that S{bi) = bi is a vertical disk in * N — N. 

Since Vo,Vk G {Ni,N2} we conclude that there exists a horizontal disk 
contained in Fix (5')n| Vo | and there exists vertical disk contained in Fix (5')n| VI- 1. 
Now, the assertion is a direct consequence of Theorem |21 I 

Proof of Theorem 1151 The assertion follows from the fact that the sets jA^il 
and |A'^2| are disjoint and we can construct an arbitrary number of sequences 
satisfying assumptions of Lemma ITHI for example 

k times 



,^ F.l ,^ F.l F,l ,^ F,l F*.-l F,-l F.-l F.-l ,^ 

Ni=^ Ni =^ . . . =^ Ni =^ Hi =^ H2 ^ H3 =^ H4 =^ N2, 



for fc > 0. 
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6 How to verify covering relations with a com- 
puter assistance. 

In this section we discuss some numerical aspects of the verification of covering 
relations. 

Let N,M be a h-sets in R" such that u{N) = u{M) = u and s{N) = 

s{M) = s and let / : |A^| R" be a continuous. In order to prove that 

f 

the covering relation N M holds, it is necessary to find the homotopy 
h : [0, 1] X A^c ^ X and a map A : R" -> R" satisfying conditions 

Since in our example fc = cm ° f ° is a diffeoniorphism we can try to 
find a homotopy between fc and its derivative computed in the center of set and 
projected onto unstable directions, i.e., we define 

A : R" 3 p ^ 7r„p/,(0)(p,0)) «E R^ 

where 7r„ : R" R" is a projection onto first u variables. We require A to be 
an isomorphism, then we have 

degiAMO, 1), 0) = sgn (det A) = ±1. 

Now we define the homotopy between fc and {p, q) {A{p),0) by 

= (l-t)/c + t(A(p),0), for e 5^(0, 1) X 57(0,1). (67) 

Obviously the homotopy H67|l satisfies conditions (|2Jl and (|SJl. We need to check 
if the homotopy H()7|l satisfies conditions H3I4|) . 
Below we describe precise algorithms. 

Definition 14 Let U C R" be a bounded set. We say that Q C 2^ is a grid of 
U if 

1. Q is a finite set 
2- UcUaegG 

3. each G E Q can be represent in a computer 

Definition 15 Let U C R 6e a bounded set. By (U)i we denote the interval 
enclosure of the set, i.e., the set (U)i is the smallest representable interval 
containing U or [—oo, oo] if there is not a representable interval containing U . 

Let f/ C R" be a bounded set. By {U)i we denote (tti ([/))/ x ■ • ■ x (7r„ ([/))/ 
where tt^ is a projection onto i-th variable. 

In the algorithms presented bellow all computations are performed in interval 
arithmetic |Mo| . 

First we discuss how we check condition 

Algorithm 1 
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function ComputeUnstableWall(5i : grid, G2 ■ grid) : bool 
var 

X, LX, Z : representable sets; 
begin 

for each Gi ^ Gi 
begin 

X Gi X 0; // e 

LX := {DMO){X))i; 

foreach G2 G G2 

begin 

Z := (/c(Gi xG2)ULX)i; 

if not TTu{Z) C M" \ S„(0, 1) return False; 

end; 

end; 

return True; 

end. 

Lemma 17 Assume N, M be a h-sets and f : |iV| R" be such that fc is 
smooth. Let Gi be a grid of dBu{0, 1) and let G2 be a grid of Bs{0, 1). If Algo- 
rithm^ is called with arguments {Gi,G2) and returns True then the homotopy 
defined in j6'7| ) satisfies condition (0)- 

Proof: Let {p,q) G . Since Gi is a grid of 9i3„(0, 1) and G2 is a grid of 
B7(0, 1) then 

Gi X ^2 := {Gi X G2 I Gi e GuG2 e ^2} 

is a grid of . Therefore {p,q) G Gi x G2 for some Gi G Gi, G2 G ^2- Since 
the Algorithm stops and returns True the condition 

MifciGi X G2)ULX),) cM"\^(0,l) 

is satisfied, which impUes that for t G [0, 1], h{t,p, q) ^ Mc- I 

Now we discuss how we verify condition Q). The main point of our approach 
is that it is enough to compute fc{dNc). 

Algorithm 2 

function ComputeBoundary(5 : grid) : bool 
var 

X : representable set; 
begin 

foreach G E G 
begin 

X := (/,(G)),; 

if not iis{X) C int_Bs(0, 1) return False; 
end; 

return True; 

end. 
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Lemma 18 Assume N, M be a h-sets, f : M" — > M" be such that fc is a 
diffeomorphism. Let Q be a grid of dNc- If Algorithmic is called with argument 
G and returns True then the homotopy defined by J6'7| ) satisfies condition 

Proof: Let x £ dNc. Since Algorithm stops and returns True we obtain that 
'^sifcix)) G 5^(0, 1). From Eq. ^ it foUows that for x G dNc and t e [0, 1] 
TTsiHt, x)) e BsiO, 1). Therefore for x 6 dNc and t e [0, 1] 

hit, x) i 5^(0, 1) X 9B,(0, 1) = Af+. (68) 

There remains to prove that Eq. H68|l is satisfied for x ^ Nc and t e [0, 1]. Since 
fc is a diffeomorphism, the Brouwer-Jordan Theorem imphes that TTs{fc{Nc)) C 
B^(0, 1). Hence, for t G [0,1] 

TTs{h(t,Nc))^Bs{0,l) 

which imphes that h{t, Nc) n M+ =0. I 
6.1 Technical data. 

The grids used in the numerical proof of Lemma 1111 and Lemma 1131 always 
consist of "boxes", i.e., products of representable intervals. The total number 
of boxes used in the proof is approximately 2.2 • 10*. The numerical proof of 
Lemma ITTI and Lemma [T^ took approximately 36 minutes on 2.4GHz processor 
under PLD Linux Distribution. 

The CH — h sources with a short description how to run the program are 
available at |W1| . 
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